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ABSTRACT 

Using a sample of more than 6000 quasars from the Sloan digital sky survey (SDSS) we 
compare the black-hole mass distributions of radio-loud and radio-quiet quasars. Based on 
the virial black-hole mass estimator the radio-loud quasars (RLQs) are found to harbour sys- 
tematically more massive black holes than radio-quiet quasars (RQQs) with very high sig- 
nificance (> 99.99%), with mean black-hole masses of < log(A/fcft/ Mq) >= 8.89 ± 0.02 
and < log(Afft?,/MQ) >= 8.69 ± 0.01 for the RLQs and RQQs respectively. Crucially, 
the new RLQ and RQQ samples have indistinguishable distributions on the redshift-optical 
luminosity plane, excluding the possibility that either parameter is responsible for the ob- 
served black-hole mass difference. Moreover, this black-hole mass difference is shown to be 
in good agreement with the optical luminosity difference observed between RLQ and RQQ 
host galaxies at low redshift (i.e. AMhost = 0.4 — 0.5 magnitudes). Within the SDSS sam- 
ples, black-hole mass is strongly correlated with both radio luminosity and the radio-loudness 
TL parameter (> 7cr significance), although the range in radio luminosity at a given black- 
hole mass is several orders of magnitude. It is therefore clear that the influence of additional 
physical parameters or evolution must also be invoked to explain the quasar radio-loudness 
dichotomy. 
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1 INTRODUCTION 

Ever since the discovery that only a small minority of optically- 
selected quasars are also luminous radio sources (e.g. Lsghz > 
10'^*WHz~^sr~^), many studies have attempted to isolate the phys- 
ical mechanism underlying this so-called quasar radio-loudness di- 
chotomy. Although previous studies found a clear bimodality in the 
radio luminosities of optically-selected quasars (eg. Kellermann et 
al. 1989; Miller, Peacock & Mead 1990), more recently the very 
existence of the radio-loudness dichotomy has been questioned 
(Lacy et al. 2001; Cirasuolo et al. 2003; although see Ivezic et 
al. 2002 for an alternative viewpoint). The principal reason behind 
this renewed interest in the radio properties of optically selected 
quasars is the ability to combine the large SDSS and 2dF opti- 
cal quasar samples with wide-area radio surveys such as the Faint 
Images of the Radio Sky at Twenty-cm (FIRST) and the NRAO 
VLA Sky Survey (NVSS). The FIRST survey (Becker, White & 
Helfand 1995) in particular has identified large numbers of so- 
called radio-intermediate quasars which have largely fiUed-in the 
apparent gap in radio luminosities between the RLQ and RQQ pop- 
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ulations (eg. Lacy et al. 2001). Consequently, the distribution of 
optically-selected quasars on the optical-radio luminosity plane is 
undoubtedly more continuous than was previously thought. Irre- 
spective of this, the fundamental question of what causes luminous 
quasars with seemingly identical optical properties to differ in their 
radio luminosities by several orders of magnitude remains unan- 
swered. 

Recent progress has been made in largely eliminating two pa- 
rameters which were originally suspected of influencing the radio- 
loudness dichotomy; host-galaxy morphology and cluster environ- 
ment. Thanks largely to the Hubble Space Telescope (HST), quasar 
host-galaxy morphologies have now been investigated out to inter- 
mediate redshifts, 0.1 < z < 0.5 (eg. Dunlop et al. 2003; Schade 
et al 2000; McLure et al. 1999; Disney et al. 1996). Drawing to- 
gether the results of these studies, it is clear that the hosts of op- 
tically luminous quasars (i.e. Mr < —24) are bulge-dominated, 
spheroidal galaxies irrespective of radio luminosity (Dunlop et al. 
2003; Schade et al. 2000). In the light of the discovery of the corre- 
lation between black-hole and bulge mass (Magorrian et al. 1998; 
Gebhardt et al. 2000; Ferrarese & Merritt 2000), this result is per- 
haps not surprising. However, it should also be remembered that in 
the optical the hosts of RLQs are consistently found to be ~ 0.5 
magnitudes brighter than their RQQ counterparts (eg. Dunlop et 



© 0000 RAS 



2 R.J. McLure and M.J Jarvis 



al. 2003), the implications of which will be discussed in Section 
3. In addition, it now appears that cluster environment does not 
greatly influence quasar radio luminosities either. Recent studies of 
the cluster environments of AGN have tended to show that quasars 
prefer to inhabit poor clusters (e.g. Abell class 0) or groups, rather 
than the cores of rich clusters (Best 2004; Sochting, Clowes & 
Campusano 2004). Where a direct comparison of matched sam- 
ples of RLQs and RQQs has been performed (eg. Wold et al. 2001; 
McLure & Dunlop 2001a) little difference has been found in their 
respective cluster environments. 

Following the inconclusive nature of the results regarding the 
Mpc and Kpc-scale environments of RLQs and RQQs, attention 
has now become focussed on possible differences in the proper- 
ties of their central engines. At the same time the robustness of 
the so-called virial black-hole mass estimator (Wandel, Peterson 
& Malkan 1999; Kaspi et al. 2000) has now become established, 
allowing estimates of AGN black-hole masses to be derived from 
a continuum luminosity and broad-line FWHM measurement only 
(i.e. Mbh oc AL'^'^FWHM^). Consequently, many studies have ex- 
ploited the virial mass estimator to investigate whether black-hole 
mass could be the hidden parameter controlling the radio-loudness 
dichotomy; with widely contrasting results. Several authors (eg. 
Laor 2000; McLure & Dunlop 2001b,2002; Dunlop et al. 2003) 
have concluded that RLQs have systematically higher black-hole 
masses than their RQQ counterparts, and have suggested the exis- 
tence of a black-hole-mass threshold (e.g. Mth — 10* Mq) above 
which the fraction of RLQs increases significantly. This separation 
of the RLQ and RQQ populations in terms of black-hole mass is 
further supported by the principal components analysis of RLQ 
and RQQ spectra by Boroson (2002). Moreover, by combining the 
FIRST Bright Quasar Survey (FBQS, Gregg et al. 1996) and the 
PG quasar survey (Green, Schmidt & Liebert 1986), Lacy et al. 
(2001) found that, given enough dynamic range, radio luminosity 
and black-hole mass were strongly correlated, albeit with substan- 
tial scatter. More recently, McLure et al. (2004) found a significant 
correlation between black-hole mass and integrated low-frequency 
(151 MHz) radio luminosity from a complete sample of 41 power- 
ful radio galaxies at « ~ 0.5. 

In contrast to the previously mentioned studies, several au- 
thors have concluded that no strong link exists between black- 
hole mass and radio luminosity in AGN (eg. Urry 2003; Oshlack 
et al. 2002; Woo & Urry 2002; Ho 2002). With particular rele- 
vance to optically luminous quasars, both Oshlack et al. (2002) and 
Woo & Urry (2002) highlight that the apparent existence of ex- 
tremely radio-loud flat-spectrum quasars (FSQs), with virial black- 
hole mass estimates in the range 10^ < M^h < 10* M©, ef- 
fectively rule out a strong correlation between black-hole mass and 
radio-loudness (although see Jarvis & McLure 2002 for a different 
interpretation of the FSQ data). 

In this paper we revisit the issue of whether black-hole mass 
and radio luminosity are related in quasars using the virial black- 
hole mass estimates calculated by McLure & Dunlop (2004) for 
objects from the SDSS quasar catalog II (Schneider et al. 2003). 
The aim of this paper is to determine whether RLQs and RQQs, 
matched in terms of redshift and optical luminosity, have distin- 
guishable black-hole mass distributions or not. Throughout the pa- 
per we adopt the spectral index convention fu oc i/^" and the fol- 
lowing cosmology: Hq = 70km s~^, JIa = 0.7, = 0.3. 



2 THE SAMPLE 

The sample adopted in this paper is drawn from the SDSS quasar 
catalog 2 (Schneider et al. 2003) which consists of 16714 quasars 
in the redshift interval 0.1 < z < 5.3 with Mi{AB) < -22. 
The black-hole mass estimates are taken from McLure & Dunlop 
(2004) who produced virial black-hole mass estimates for more 
than 12000 of the SDSS quasar catalog 2 (SQC2) objects with red- 
shifts z < 2.1, based on both the H/9 and Mgll emission lines. 
Full details of the modelling of the SDSS quasar spectra and the 
calibrations of the virial mass estimator can be found in McLure & 
Dunlop (2004) and McLure & Jarvis (2002). The final samples of 
RLQs and RQQs were selected from the full SQC2 via appUcation 
of the following criteria: 

(i) All objects which fell outside the SDSS/FIRST overlap re- 
gion were excluded to ensure consistent radio data for the final 
sample. 

(ii) All objects with i{AB) > 19.1 or which did not satisfy the 
SDSS quasar colour selection algorithm for follow-up spectroscopy 
were also excluded. This criterion ensures that the final samples 
are not biased by the inclusion of pre-selected FIRST sources or 
serendipitous targets (Schneider et al. 2003). 

The remaining objects were cross-correlated with the McLure 
& Dunlop (2004) black-hole mass estimates which cover > 90% 
of the SDQ2 objects with z < 2.1. This final list of quasars was 
then separated into the RLQ and RQQ samples based on the radio- 
loudness TZ parameter (TZ = fzoHz/ fs, where fsoHz ^ and /s 
are the observed flux densities at 5GHz and 4400A respectively). 
Following the standard convention we classify every object with 
TZ > 10 as radio-loud, and everything with TZ < 10 as radio- 
quiet. An upper limit in TZ was calculated for quasars undetected 
by FIRST using the nominal FIRST object-detection threshold of 
ImJy. To investigate the possibility that the quasar radio luminosi- 
ties are underestimated, due to the FIRST survey resolving out ex- 
tended flux (e.g. Blundell 2003), we cross-correlated our samples 
with the lower resolution NVSS survey. The results of this showed 
no significant differences in 1.4GHz flux densities, which suggests 
that the radio emission of these optically-selected quasars is fairly 
compact. 

The final RQQ and RLQ samples produced by our selection 
process comprise 6099 and 436 objects respectively. We note here 
that adopting TZ > 10 as our radio-loudness threshold is conser- 
vative, in the sense that 30% of the final RLQ sample have 5GHz 
radio luminosities less than the alternative radio-luminosity based 
threshold of lO^^WHz'^sr"^ (Miller, Peacock & Mead 1990). The 
results of adopting the strict radio luminosity threshold are dis- 
cussed further in Section 3. 



2.1 Matched samples 

In order to investigate any differences in the black-hole mass distri- 
butions of the RLQs and RQQs the samples should be free from 
biases which could affect the virial black-hole mass estimator. 
Given the role of optical luminosity in the virial mass estimator 
(Mth XL'^^'FWHM^), this essentially requires that the two sam- 
ples are indistinguishable in terms of their optical luminosity and 
redshift distributions. 

^ The 5GHz flux-density is extrapolated from 1 .4GHz using the mean spec- 
tral index of FBQS quasars; a = 0.5 (White et al. 2000) 
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Redshift z Log (M^i^/Mj^) 

Figure 1. (a) Bolometric luminosity versus redshift for the RQQ (open grey circles) and RLQ (filled black circles) samples. (b) Bolometric luminosity versus 
virial black-hole mass estimate for the same samples, with those RLQs with Lsghz < 10^'^WHz^^sr^^ plotted as open black circles. 



In Fig^ we show the RLQ and RQQ samples on the Li,oi — z 
plane where, following McLure & Dunlop (2004), the quasar bolo- 
metric luminosities are estimated from either the 5100A or 3000A 
luminosities as Lboi ~ 9.8AI/5ioo or Lboi = 5.9AI/3ooo, depend- 
ing on their redshift (McLure & Dunlop 2004 demonstrated that 
the ratio of these two bolometric luminosity estimators is unity for 
~ 1200 quasars where both could be measured). The application of 
the one-dimensional Kolmogorov-Smirnov (KS) test finds that both 
the Lboi and redshift distributions of the RLQ and RQQ samples 
are indistinguishable; p = 0.16 andp = 0.50 respectively. In addi- 
tion, an application of the two-dimensional KS test to the Lboi — z 
distribution shown in Fig^ confirms that the two-dimensional dis- 
tributions of the RLQ and RQQ samples are also indistinguishable; 
p = 0.21. Consequently, any black-hole mass difference between 
the two quasar samples should not be a product of significant biases 
in their redshift or optical luminosity distributions. 



3 BLACK-HOLE MASS DISTRIBUTIONS 

The distributions of the RQQ and RLQ samples on the Lboi — Mbh 
plane are shown in FigQ). The application of the two-dimensional 
KS test returns a probability of only p = 1.2 x 10"^'^ that both sam- 
ples are drawn from the same parent population. The mean black- 
hole masses for the RQQ and RLQ samples are as follows: 

RLQ: < log(M6h/M0) >= 8.89 ±0.02 (8.93) 

RQQ : < log{Mbh/ Mo) 8.69 ± 0.01 (8.74) 

where the quoted uncertainties are the standard errors and the val- 
ues in parenthesis are the sample medians. Although the optical 
luminosity distributions of the RLQ and RQQ samples are indis- 
tinguishable, the mean optical luminosity of the RLQs is 0.06 dex 
brighter than that of the RQQs. Given that the virial mass estima- 
tor adopted here has a luminosity dependence of AL^ '' (McLure 
& Dunlop 2004), this will account for 0.04 dex of the black-hole 
mass difference between the samples. Consequently, accounting for 
this slight luminosity bias, we conclude that RLQs typically har- 
bour black holes 0.16 dex (45%) more massive than their radio- 
quiet counterparts. Moreover, by calculating the mean RLQ— RQQ 
black-hole mass difference in narrow luminosity bins (0.3 dex 
width) it was confirmed that the black-hole mass difference is con- 



sistent with a constant 0.16 dex off-set over the full luminosity 
range. 

3.1 Line-width difference 

The measured difference in mean black-hole mass between the 
RLQ and RQQ samples is predominantly due to fact that the 
FWHM of the low-ionization H/9 and Mgll emission lines of the 
RLQs are ~ 20% larger than in the RQQs of the same optical lu- 
minosity. The conclusion that RLQs exhibit broader low-ionization 
emission lines than RQQs is in good agreement with many previ- 
ous studies (eg. Boroson & Green 1992; Corbin 1997; McLure & 
Dunlop 2002; Sulentic et al. 2003). However, although the velocity- 
width difference between RLQs and RQQs has been noted pre- 
viously, because of potential selection effects, attributing this di- 
rectly to a difference in black-hole mass has been problematic. For 
example, given that the FWHM of low-ionization emission lines 
in RLQs are known to be influenced by viewing angle (Wills & 
Browne 1986; Brotherton et al. 1996), the velocity- width differ- 
ence between RLQs and RQQs could potentially be an orientation 
effect produced by comparing radio and optically selected samples. 
However, as discussed in Section 2, here both samples have been 
selected consistently by the SDSS ugri colour algorithm, with tar- 
geted FIRST sources deliberately excluded. Combined with the ex- 
cellent sample matching in terms of optical luminosity and red- 
shift, the influence of potential orientation effects should therefore 
be minimized. Consequently, the simplest interpretation of these 
new results is that at a fixed optical luminosity RLQs do, on aver- 
age, harbour more massive black-holes than RQQs, accreting at a 
commensurately lower fraction of their Eddington limit. 

3.2 Comparison with host-galaxy results 

Irrespective of the virial black-hole mass estimates presented here, 
it could reasonably be argued that the correlation between bulge 
and black-hole mass (Mbh — Mbuige) and the well known K — z 
relation for powerful radio galaxies (Lilly & Longair 1984) imme- 
diately imply that RLQs should harbour black holes with masses 
>10** M0. The recent Willott et al. (2003) study showed that the 
radio-galaxy K — z relation is consistent with the passive evolution 
of a ~ 3L* elliptical formed at high redshift {z > 5). Therefore, 
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Figure 2. (a) Radio-loudness parameter TZ. versus black-hole mass, (b) Absolute 5GHz radio luminosity versus black-hole mass. Symbols as Fig 1. RQQs 
undetected by FIRST are shown as upper limits. A black-hole mass of 10* Mq is highlighted by a solid vertical line in both figures. 



if the predictions of orientation-based unification are correct, the 
host-galaxies of RLQs of comparable radio luminosity should have 
similar AT— band luminosities . 

Notably, using the McLure & Dunlop (2002) fit to the Mbh - 
Mbuige relation, the mean black-hole mass of the RLQ sample cor- 
responds to a if— band bulge luminosity of 3.8 ± 0.5L* (assuming 
a 2 = colour ofR-K ^ 2.7; Bruzual & Chariot 2003), in good 
agreement with the Willott et a\. K — z relation. Moreover, it is also 
interesting to compare the RLQ black-hole mass estimates derived 
here with those derived via host-galaxy luminosities for a complete 
sample of 41 2 ~ 0.5 radio galaxies by McLure et al. (2004). Based 
on HST imaging of radio galaxies spanning the same range of ra- 
dio luminosities as the current RLQ sample, the mean black-hole 
mass was determined to be < log(M6h/MQ) >= 8.87 ± 0.04, 
in excellent agreement with the mean black-hole mass of the RLQ 
sample derived here. 

In addition to comparisons with radio galaxy studies, it is also 
possible to compare the black-hole mass difference between the 
RLQ and RQQ samples determined here with imaging studies of 
quasar host-galaxies. The most accurate quasar host-galaxy studies 
are those based on high-resolution HST imaging data (eg. Dunlop 
et al. 2003; Schade et al. 2000; Disney et al. 1995). Among the 
various HST-based host galaxy studies, the Dunlop et al. (2003) 
study has the advantage of featuring RLQ and RQQ samples 
which have well matched redshift - optical luminosity distributions. 
Based on the results of this i?— band imaging study, the hosts of 
RLQs are typically ~ 0.4 mags brighter than their RQQ coun- 
terparts. Combining this with the McLure & Dunlop (2002) fit to 
the Mbh — Mbuige relation implies that RLQ and RQQ black-hole 
masses should differ by ~ 0.2 dex, in good agreement with the 
virial-based results derived here. 



4 THE CORRELATION BETWEEN BLACK-HOLE MASS 
AND RADIO LOUDNESS 

In Fig|2|we plot both the radio-loudness TZ parameter and abso- 
lute radio luminosity versus black-hole mass for the RLQ and RQQ 
samples. In order to assess the significance of the apparent corre- 
lations in the presence of upper limits we have applied the gener- 
alized Kendall's r test in the ASURV package (Isobe, Feigelson 



& Nelson 1986) . The Kendall's r test detects a highly significant 
correlation in both cases, returning values of r = 0.0497 (7. 60-) 
and T = 0.0697 (11.4cr) for the 7^ - Mbh and Lsghz - Mbh 
correlations respectively. 

The detection of a significant correlation between black-hole 
mass and radio-loudness is in good agreement with the Lacy et 
al. (2001) study of the combined FBQS-l-PG sample. Furthermore, 
it can be seen from Fig [2^ that the RLQs are virtually exclu- 
sively confined to Mbh > 10* Mq, as previously found by Laor 
(2000), Lacy et al. (2001), Boroson (2002) and McLure & Dunlop 
(2002). Qualitatively, the fraction of RQQs with Mth < 10** Mq 
is 9.2% ± 0.4, whereas the equivalent fraction of RLQs is only 
3.7% ± 1.0. In fact, if the definition of radio-loud is restricted to 
only those quasars with Lsghz > lO^^WHz-^sr-^ only two ob- 
jects have Mbh < 10* Mq (< 1%). This result, combined with 
the radio galaxy K ~ z relation (e.g. Jarvis et al. 2001; Willott et 
al. 2003), strengthens the claim that regardless of the exact nature 
of the black-hole mass - radio luminosity relationship in luminous 
quasars, sample selection based on high radio luminosity is effec- 
tive in isolating the largest black holes, and presumably host galax- 
ies, at all epochs (eg. McLure 2003). 

4.1 Flat-spectrum quasars 

In contrast to the results found here, we note that Woo & Urry 
(2002) found no indication of a 7^ — Mbh correlation in their study 
of a heterogeneous sample of 747 quasars in the redshift interval 
< z < 2.5. The simple reason for this apparent discrepancy is 
that within their sample Woo & Urry include the 39 flat-spectrum 
quasars (FSQs) studied by Oshlack et al. (2002). Jarvis & McLure 
(2002) have examined the likely effects of orientation on both the 
radio luminosities and black-hole mass estimates for the Oshlack et 
al. (2002) sample. Although estimates of the intrinsic unbeamed TZ 
parameters are hard to accurately determine, Jarvis & McLure con- 
cluded that the radio luminosities and black-hole mass estimates 
of Oshlack et al. (2002) are likely to be over/under-estimated by 
factors of ~ 100 and ~ 4 respectively. 

Consequently, in all likelihood the Oshlack et al. FSQs are 
not inconsistent with the Lsghz — Mbh correlation seen in Fig 
with all the FSQs with intrinsic radio luminosities Lsghz > 
lO^^'WHz^^sr"^ likely to have Mi,h> 10* Mq. Notably, if the 
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FSQs are removed from the Woo & Urry (2002) sample, then their 
resuks are entirely consistent with those derived here. 



5 CONCLUSIONS 

In this paper we have compared the black-hole mass distributions 
of RQQs and RLQs using large samples which have indistinguish- 
able distributions in the redshift-optical luminosity plane. The re- 
sults of this comparison demonstrate that at a fixed optical lumi- 
nosity RLQs harbour black-hole masses which are typically 0.16 
dex (45%) more massive than their radio-quiet counterparts. In 
agreement with previous studies, this black-hole mass difference 
is shown to be largely produced by RLQs displaying broader low- 
ionization emission lines than comparable RQQs. Moreover, using 
the combined RQQ-l-RLQ sample we have shown that black-hole 
mass is strongly correlated with both absolute radio luminosity 
and the radio-loudness parameter TZ. Furthermore, the black-hole 
masses of genuine RLQs are found to be virtually exclusively con- 
fined to Mhh > 10* Mq, in good agreement with the conclusions 
of previous studies (Laor 2000; Boroson 2002, McLure & Dunlop 
2002). 

Combining the virial black-hole mass estimates with the local 
Mbh — Mbuige relation suggests that RLQ host galaxies are fully 
consistent with the radio galaxy K — z relation, in agreement with 
AGN orientation-based unification. It is therefore clear that sam- 
ple selection based on high radio luminosity is an effective method 
for cleanly isolating the largest black-holes, and presumably host 
galaxies, at all epochs. However, the new results presented here also 
confirm that no tight correlation exists between radio luminosity 
and black-hole mass, and that at a fixed black-hole mass the range 
in quasar radio luminosities covers several orders of magnitude. 
Consequently, there is clearly a large overlap between the RLQ and 
RQQ populations, with many RQQs being indistinguishable from 
RLQs in terms of black-hole mass, accretion rate, orientation, host 
galaxy properties and cluster environment. In conclusion, it there- 
fore appears inescapable that the influence of additional physical 
parameters such as black-hole spin and/or evolutionary effects must 
contribute to the radio-loudness dichotomy. 
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